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Abstract 

We study the Anderson model of localization with anisotropic hopping in 
three dimensions for weakly coupled chains and weakly coupled planes. The 
eigenstates of the Hamiltonian, as computed by Lanczos diagonalization for 
systems of sizes up to 48^, show multifractal behavior at the metal- insulator 
transition even for strong anisotropy. The critical disorder strength Wc de- 
termined from the system size dependence of the singularity spectra is in 
a reasonable agreement with a recent study using transfer matrix methods. 
But the respective spectrum at Wc deviates from the "characteristic spec- 
trum" determined for the isotropic system. This indicates a quantitative 
difference of the multifractal properties of states of the anisotropic as com- 
pared to the isotropic system. Further, we calculate the Kubo conductivity 
for given anisotropics by exact diagonalization. Already for small system sizes 
of only 12^ sites we observe a rapidly decreasing conductivity in the directions 

with reduced hopping if the coupling becomes weaker. 
71.30.+h, 72.15.Rn 
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I. INTRODUCTION 



It is well known that the three dimensional (3D) isotropic Anderson model exhibits 
a metal- insulator transition (MIT): Increasing the disorder of the random potential site 
energies causes the wave functions to localize.§ There exists a mobility edge in the energy- 
disorder diagram which separates extended from localized eigenstates. In order to determine 
these critical disorders Wc{E) accurately, the transfer- matrix method (TMM) together with 
the one-parameter finite-size scaling hypothesis applied to quasi-lD bars has been used with 
much success in the past.i^ Recently, the anisotropic Anderson model has received much 
attention in connection with the anisotropic transport properties of the high Tc cuprates 
and a possible contradiction to the scaling theory of localization was mentionedi supported 
by a diagrammatic analysis^ However, recent TMM studiesi^ show that the one-parameter 
scaling theory is still valid and further that an MIT exists even for strong hopping anisotropy 
7. The values of the critical disorder in the band center were found to follow a power law 
Wc oc (1 — 7)^ independent of the orientation of the quasi- ID bar. P was argued to be 
independent of the strength of the anisotropy. 

Here, we shall study the problem of Anderson localization by a different method: we 
focus our attention directly on the eigenf unctions of the Hamiltonian. In an infinite system 
the wave functions are expected to be localized on the insulating side and extended on the 
metallic side even arbitrarily close to the transition. As first suggested by Aoki§ the fractal 
nature of the critical eigenstates can connect these discrepant characteristics. Indeed, large 
fluctuations of the wave functions have been observed numerically which dominate — at 
least at small length scales — the character of the states and invalidate the simple notions 
of exponentially localized or homogeneously extended states. Approaching the transition 
these fluctuations increase and at the critical disorder they are expected to occur on all 
length scales. It has been shown0 that such wave functions are multifractal entities. To 
characterize the eigenstates of the isotropic Anderson model the singularity spectrum f{a) 
has been used successfully.0 A characteristic spectrum /c(a) was shown to determine the 
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mobility edge independent of the microscopic details of the sample. l3 Further, around its 
maximum, /c(a) agrees well with an analytical result of Wegneilll based on a nonlinear 
a model calculation. Near the critical disorder Wc-, characteristic changes of /(«) were 
observed when the system size was incteased.Q These distinguish the localized and the 
extended character of the states and therefore allow us to determine the transition directly 
from multifractal properties of eigenstates. 

It is our aim in the present work to use and extend these concepts for the case of 
anisotropic hopping. In Sec. J| we introduce our notation and define the anisotropics of 
weakly coupled planes and weakly coupled chains. We next recall the concepts and methods 
of the multifractal analysis employed in the sequel. Using the hypothesis of a characteristic 
singularity spectrum, we estimate the critical disorders in Sec. [IV B| . To check the validity 
of the hypothesis we analyze the system size dependence of the multifractal properties and 
compare our results with TMM dataSi For completeness, we also study the conductivity of 
small samples of anisotropic systems in Sec. |V|. We discuss our results in Sec. [V^. 

II. THE ANDERSON MODEL WITH ANISOTROPIC HOPPING 

The Anderson Hamiltonian is given as0 

i ijtj 

Here, the sites i = {x, y, z) form a regular cubic lattice of size A^^ and the potential site 
energies are as usual taken to be randomly distributed in the interval [—W/2,+W/2]. 
The transfer integrals tij are restricted to nearest neighbors and depend only on the spatial 
direction, so t^j can either be t^^, ty or t^. We set the energy scale by normalizing the largest 
tij to 1. 

Following Ref. P, we study two possibilities of anisotropic transport: (i) weakly coupled 
planes with tx = ty = 1^ tz = 1 — "i and (ii) weakly coupled chains with t^ = ty = 1 — j, 
tz = 1. Here the parameter 7 G [0, 1] describes the strength of the anisotropy. Hence, for 
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7 = we recover the isotropic 3D case and 7 = 1 corresponds to independent planes or 
A^^ independent chains. The direction with normal (reduced) transfer integral is called the 
parallel (perpendicular) direction. 

The Lanczos algorithm,lii which is well suited for the diagonalization of sparse matrices, 
allows us to solve the eigenvalue equation H^{E) = E'^{E) for system sizes up to = 72, 
yielding eigenvalue/eigenvector pairs in a requested energy range. We use state-of-the-art 
workstations and a parallel computer with 128 PowerPC processors. It takes about 11 
hours to diagonalize the Anderson Hamiltonian with A^ = 48 on the parallel machine using 

8 processors. The workstations need about 35 hours for the same calculation. Since we 
also have to perform a statistical averaging over different disorder configurations, we have 
restricted the systematic investigations to sizes up to A^ = 48. In order to allow a direct 
comparison with the results of Refs. 0-^, we restrict our study to the states in the center of 
the band such that E = Numerically this is the worst case because of the high density 
of states there which requires a very large tridiagonal matrix in the Lanczos algorithm to 
determine the eigenvalues. 



III. MULTIFRACTAL ANALYSIS 

Fractal measures are widely used in physics to characterize objects such as percolat- 
ing clusters, random walks, and random surfaces, mil The common geometric feature of 
such point sets is the self- similarity: Parts of the set are similar to the whole, at least in 
a statistical sense. However, for fluctuating physical quantities such as the probability am- 
plitude of an eigenstate \E'(i?) of the Anderson model, the appropriate concept is given by 
the multifractal measures: If the mentioned fluctuations are statistically the same on every 
length scale, i.e., if all the moments of the investigated quantity are self similar, the object 
is (statistically) self-afiine and is called a multifractal. 

A characteristic property of multifractals is their singularity spectrum /(a).0 Let us 
briefly describe an algorithm to determine this quantity, based on the standard box counting 
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procedure: We consider a volume L in our D dimensional space which contains the support 
of the physical variable, i.e., all points where the variable is defined. We cover it with a 
number of "boxes" of linear size r = 6L. The actual shape of the boxes is not important, 
they may be spheres as well. Next, we determine the contents yUi(5) of each box i by summing 
or integrating the investigated quantity over the part of the support inside the box. For a 
self-affine object one finds a power-law dependence fii{S) oc 6°'' in the limit 5 — > 0. The so- 
defined singularity strength is assigned to each point of the support. The subset Sa which 
contains all points with the same value of a is a fractal with fractal dimension /(a) defined 
by K{a,6) oc 6~-^^°'\ Here, K{a,6) is the number of boxes which cover Sa- A multifractal 
object consists of a (infinite) number of subsets Sa with different fractal dimensions. 

In the present work we shall use an equivalent but numerically more convenient algo- 
rithm to compute the singularity spectrum. Our physical quantity is again the probability 
amplitude of eigenstates. Considering the normalized gth moments of the box probability 
6) = J2k lA{^) it possible to find^l a parametric expression of /(a) such that 

= lim5^oEi/^i(g,'5)lnAii(l,'^)/ln5 , 

(2) 

/(g) = lim^^o Ei /^i(g, 5)ln/ii(g, 5)/\\i5 . 

We plot the sums in @j versus ln5 and observe multifractal behavior if and only if the data 
may be well fitted by straight lines. The slope from the linear regression procedure used in 
the fit gives / and a. Note, that a check of the linearity is important, since the numerical 
procedure gives an f{a) curve for nearly every distribution of the physical variable, but 
without the mentioned linearity it does not indicate multifractality. 

In general, /(a) is a nonnegative, convex function with < amin < « < Oimax < oo. The 
maximum of /(a) at a{q = 0) = cto equals the dimension of the support, i.e., the fractal 
dimension -D/ of the subset of points where the investigated quantity is not zero. For our 
wave functions Dj = D = 3 because they are nowhere exactly zero. The whole /(«) curve 
is below the bisector f{a) = a except at a{q = 1) = ai where both curves touch. For 
q = 1 the relation /(cti) = cti is fulfilled, ai equals the entropy dimension or information 
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dimension and one can show that the corresponding set S'qj contains the entire measure.tS 
There are two hmits which will be important for the later interpretation of our results. 
Consider a L'-dimensional support, (i) A uniform distribution is represented by the single 
point f{a = D) = D in the singularity spectrum, because fii{6) oc {L6)^ for every point of 
the support, (ii) With increasing localization the spectrum becomes wider and an extremely 
localized distribution with measure 1 at one point and elsewhere has a spectrum which 
consists of two points only: /(a = 0) = and /(a = oo) = D. This is because the box 
around the maximum has contents 1 for each 6, so a is for this single point while all 
other points have fii (x 6°° = 0. In Fig. |I] we show two typical singularity spectra of 3D 
wave functions corresponding to a localized and an extended wave function. The tendency 
towards the 2 limiting cases can be seen for these two examples already: The extended state 
has a narrow /(a) curve close to /(3) =3 while the localized wave function is represented 
by a very wide spectrum with larger ao and smaller ai. 



IV. CALCULATION OF CRITICAL DISORDERS Wc{j) 

A. Existence of multifractal eigenstates 

As has been shown in Refs. by the TMM, the anisotropic Anderson model still 
exhibits a MIT for all 7 > in the band center E = and, by the general arguments given 
above, we expect the wave functions at the transition point to be multifractals just as in 
the isotropic case. As a check we have computed various eigenstates close to the proposed!'! 
critical disorders Wc for system sizes up to = 48. In Fig. ^ we show the data for the 
linear regression of a typical state with W ~ Wc- We find even for very strong anisotropies 
7 = 0.99 that the sums in Eq. (0), plotted versus InS are linear. Therefore, we do find 
multifractal behavior of the wave functions close to the critical disorder for the anisotropic 
Anderson model. 

Every singularity spectrum is characteristic only for the particular configuration of the 
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site energies. But for a given set of parameters {W, E, 7} the different /(a) curves fluctuate 
around one singularity spectrum. In order to suppress these statistical fluctuations we 
average the spectra obtained from 3 to 8 eigenstates close to E = for 12 realizations of the 
random site energies. The averaged spectrum is thus characteristic for the set of parameters 
{W,E,'y} and will be used in the next sections to compute the critical disorder Wc as a 
function of the anisotropy 7. 

B. Estimation of Wc from comparison with the characteristic spectrum 

In the isotropic case a characteristic singularity spectrum fc{a) was found previouslyil 
at all points of the mobility edge independent of the microscopic details of the system such 
as the probability distribution of the site energies. The region close to the maximum of 
/c(q;) is described well by an analytical result of WegneJii from the 2 + s expansion of the 
non- linear a model, i.e., 



As a hypothesis we shall now assume that this characteristic spectrum determines the tran- 
sition even in the case of anisotropic hopping. This hypothesis is certainly valid in the limit 
7-^0 but needs further support for larger anisotropics. 

We find that for each anisotropy 7 there exists a corresponding Wf^ such that the eigen- 
states are characterized by fc- Identifying Wc = Wf^ gives us an estimate for the 7 depen- 
dence of the critical disorder. Note that since the singularity spectrum should be independent 
of the system size at the transition point, it is sufficient to investigate small systems. We 
have used systems with N = 24 for the results presented in this section. 



Assuming the validity of fc we find a crossover between two power laws in the 7 depen- 
dence of the critical disorder: Wc = 55(1 — 7)°-^^ for 7 > 0.9 and Wc = 16.8 (1 — 7)°-^^ for 
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1. Weakly coupled planes 



7 < 0.9 as can be seen in Fig. ^. This does not agree with the results of Ref. ^, where 
13 = 0.25 has been calculated within the self-consistent theory of localization and where the 
single power law Wc = 15.4(1 — 7)°-^^ has been deduced from the TMM data. 



2. Weakly coupled chains 

In Fig. ^ the results for Wc{'~f) for weakly coupled chains are shown. Using we find 
Wc = 17.6(1 - 7)°-^^ which is very similar to the TMM dat J Wc = 16.19(1 - -ff-<^^\ The 
difference becomes significant only for very large 7 ^ 0.9. The exponent (3 = 0.611 was 
obtained from a fit of the TMM data over the whole 7 range. For large 7 the authors of 
Ref. I get P = 0.5. This is consistent with the result of Ref. ||. 



C. Estimation of Wc from the system-size dependence 

We have shown in the last section, that the assumption of the characteristic fc leads to 
large differences in the estimates of Wc between the TMM results&i and our results based 
on the multifractal analysis. Thus we will now use a more direct method to estimate Wc{'~f) 
from the multifractal properties of the eigenstates. From the isotropic case it is knownlll 
that multifractal behavior is found not only directly at the critical disorder Wc but also 
close to the transition. The reason is the finite sample size which is much smaller than the 
characteristic length scales of the states close to Wc- In this range the exponential decay 
or uniformly extended character of the wave function is masked by large fluctuations and it 
is not obvious to which side of the MIT a given state belongs. Due to the relatively small 
sample size the system is very sensitive to its boundary. Correspondingly, a characteristic 
change in the singularity spectrum is observed when the system size is increased. This change 
can be evaluated to distinguish the localized or extended character of the wave function. For 
an extended state the /(a) curve becomes narrower and the maximum position is shifted 
towards smaller values of a, approaching the value 3. The opposite behavior is found for 



a localized state. Thus the spectra tend towards the extreme cases discussed in Sec. |T1 



Indeed we expect these limiting cases, namely /(3) = 3 for the metallic side and /(O) = 
and /(oo) = 3 for the insulating side, to be reached for infinitely large system size for any 
disorder except Wc- Only directly at the transition the wave functions are multifractal, 
the fluctuations are the same on all length scales, and /(a) is independent of the system 
size. This makes is feasible to determine the critical disorder by analyzing the system-size 
dependence of the singularity spectra.0 



1. Weakly coupled planes 

We show in Fig. ^ an example of weakly coupled planes with 7 = 0.8. The above 
described different behaviors of the spectra can be seen. For W = 8, a larger system size 
results in a narrower /(a) curve which is characteristic for extended states. On the other 
hand, the increase in the system size for W = 12 yields a widening of the spectrum indicating 
localized states. The singularity spectrum for = 10 is least effected by the change of the 
system size and we thus conclude a critical disorder Wdj = 0.8) = 10± 1. Considering the 
error bars, the /(a) curve for Wc equals the characteristic spectrum fc of the isotropic case. 
For moderate anisotropics 7 ^ 0.8 this confirms the hypothesis of a characteristic fc- 

Visual observation of the system-size dependence of the /(a) curves is not well suited 
for a systematic search for the transition. A better method is to focus attention to special 
points of the spectra such as the position ao of the maximum and the information dimension 
Di = ai. An increase of the system size causes a decreasing ao and a increasing ai for 
extended states and the opposite tendency for localized states0ill as described in Sec. |T]. 



A constant behavior of and ai versus system size indicates Wc- Following Ref. |TJ we 
have parameterized the system-size dependence by l/ln(A^) which has been found to give a 
nearly linear behavior of ao and ai thus distinguishing their tendencies more clearly.0iil 
In Fig. ^ we find a constant behavior at the same value of the disorder for both quantities 
and we conclude Wd'j = 0.96) = 8.0 ± 0.5. 



For very weakly coupled planes we get significantly larger values of Wc than in Sec. [IV B 
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The new values are close to, but slightly larger than the TMM data! as can be seen in 
Fig. ^. Our data follow Wc = 16.3(1 — 7)'^-^^ which confirms the exponent (3 = 0.25 derived 
analytically.^ We therefore conclude that fc is no longer characteristic for the eigenstates 
at the MIT for weakly coupled planes with 7 ^ 0.8. In our present analysis we find wider 
singularity spectra which is a sign of a tendency towards localization. An eigenstate at the 
transition for very strong anisotropy 7 = 0.99 is shown in Fig. ^ The probability amplitude 
is concentrated to a few planes perpendicular to the direction with reduced transfer. This 
coincides with the observation that the localization length is smaller by a factor 1 — 7 in 
the perpendicular direction compared with the parallel one.i In the mentioned planes the 
wave function looks like a fractal with holes and islands of different sizes, very similar to 
the critical eigenstates of the isotropic system. It may well be that the cubic boxes used in 
the box-counting procedure for the multifractal analysis cannot appropriately measure this 
fractal, because most box sizes exceed the number of planes on which the wave functions 
are concentrated. Therefore it is possible that the deviations of Wf^ from Wc in Fig. ^ are 
an artefact of the analysis. 



2. Weakly coupled chains 

The results for the Wc{'y) dependence of weakly coupled chains are shown in Fig. 
They are in reasonable agreement with the TMM data,! although we cannot reproduce the 
exponenti f3 = 0.5. The differences between Wc and Wf^ are not as large as in the other 
case and the multifractal properties of the critical states are therefore similar to those of the 
isotropic system. 



V. CONDUCTIVITY IN SMALL ANISOTROPIC SYSTEMS 

The transport properties are determined by the localization properties of the states. 
At T = localized states cannot contribute to charge transfer and we have insulating 
behavior. On the other hand, extended states yield metallic behavior. The Kubo formula 
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following from Fermi's golden rule provides a connection of the electrical conductivity and 
the electronic states \n). 

Let us consider an electrical AC field with frequency uj = hE in x direction on an sample 
with volume A^^. We assume a half-filled band at T = such that all states with E < are 
occupied while all others are unoccupied. A configuration average because of the random 
site energies is denoted by ( )^. Neglecting prefactors the real part a of the conductivity is 
given byii 



In order to compute this quantity it is necessary to know all eigenvalues and all eigen- 
states. We use the standard Householder algorithm to diagonalize the Hamiltonian (p. 
We average over 90 to 150 configurations to suppress the large statistical fiuctuations of 
a. This limits the treatable linear system size to A^ = 12. As a consequence we encounter 
strong finite-size effects for W < A. The small number of eigenenergies in the ordered 
limit is not smeared out sufficiently by the disorder to yield a smooth density of states 
p{E) = {J2n ^{E - En))(j as shown in Fig. ^ for = 1 and 7 = 0.9. We also note that 
the density of states for larger disorder values and 7 > 0.9 agrees with that of the respective 
ID or 2D system within the uncertainty due to fiuctuations. However, the transport behav- 
ior of the states is completely different: In ID and 2D there is no MIT {Wc = 0) while 3D 
systems exhibit an MIT even for very strong anisotropy as shown in Sec. 

Because the characteristic length scales of the wave functions exceed the system size it is 
a priori not clear whether the localization behavior of the states has any measurable influence 
on the computed conductivity. Suppose there is no such influence, then the matrix elements 
{n\x\n') are all be equal and the conductivity is given by the i?- weighted joint density of 
states au{E) = {J2nn' Enn'S{E — Enn'))c- Wc Compare a and o"u for weakly coupled planes 
with, e.g. W = 1 and 7 = 0.9 in Fig. |^. The peak structure for this small disorder is a 
again flnite size effect reflecting the peaks of p{E). The positions of the minima of 
are the same as expected from the joint density of states but the minima are much more 




(4) 
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pronounced. The reason for this behavior is the strong locahzation of the states for energies 
with low p{E) similar to the localization in the band tails; the latter causes the decrease of a 
at higher energies. Thus despite the small system size the conductivity is highly influenced 
by localization effects. 

In Fig. |T^ we present the conductivity computed for the two nonequivalent directions: 
parallel (a), (c) and perpendicular (b), (d) to the planes and chains, respectively. For 
strong anisotropy the wave functions are concentrated to a few chains or planes as shown in 
Fig. ^ Consequently, the conductivity is drastically reduced in the perpendicular direction. 
For 7 = 0.9 the maximum of a is reached at small energies because the most extended 
eigenstates appear in the band center. This causes the (small) peak for W = 1. For strong 
disorder W = 15 all eigenstates are strongly localized and the perpendicular conductivity is 
nearly negligible. For the parallel conductivity we find an increase if the anisotropy becomes 
stronger. Here, the transport is not handicapped by the anisotropic localization of the wave 
functions. The increase is relatively small for the planes and considerable for the chains. A 
good argument to explain this difference is the form of the density of states which yields a 
higher amount of possibly transitions for the energies around the position of the maximum of 
cr{E) for the second. The parallel conductivity at = 15 is relatively small but considerable 
in a large energy range which reflects the disorder widened energy band. We note that the 
conductivities for a very strong anisotropy 7 = 0.99 are nearly equal to those of 7 = 0.9 in 
the parallel direction and again negligible in the perpendicular direction. 

VI. CONCLUSIONS 

In the present work we have studied the localization behavior of eigenfunctions and 
transport properties of the Anderson model with anisotropic hopping. As expected from the 
general argument for the fractal nature of wave functions at the metal-insulator transition, 
multifractal eigenstates were found even for strong anisotropy. The multifractal description 
holds not only directly at the transition but also close to it due to the small system sizes 



12 



considered. As a first estimate for tlie critical disorder Wc we determined tliat Wf^ wliere 
tlie states sliow tlie cliaracteristic singularity spectrum fda) which indicates the MIT in the 
isotropic case. But especially for weakly coupled planes the computed anisotropy depen- 
dence of the critical disorder differs remarkably from the TMM results.i We also analyzed 
the system-size dependences of the singularity spectra to determine the MIT. The observed 
Wc{'j) agree reasonably well with the TMM data.i'i Therefore we conclude that the "char- 
acteristic spectrum" is no longer valid if the anisotropy becomes strong. This is surprising 
because fc was independent of the microscopic details of the isotropic system for the 3D 
case. The spectrum at Wc for weakly coupled planes is wider than fc- This coincides with 
the observed concentration of the probability amplitude to only a few planes perpendicular 
to the direction with reduced hopping for large anisotropy. 

We have also studied the AC conductivity of small anisotropic samples using Kubo's 
formula. In this case the treatable system sizes are very small because all eigenvalues and 
eigenvectors of the Hamiltonian are needed. Nevertheless we observe a rapidly decreas- 
ing conductivity in the direction with smaller hopping integral if the anisotropy becomes 
stronger. This is a pure localization effect. For the used small system size = 12 it 
is surprising that this can be observed, because the characteristic length scales are much 
larger for nearly all of the states. Another interesting fact is that the density of states for 
an anisotropy 7 = 0.9 is already nearly identically with that of the corresponding lower- 
dimensional system. 
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FIG. 1. Singularity spectra of a localized {W = 1) and an extended {W = 25) state of an 
isotropic system with A'' = 48. The circles (o) mark /(ao) and the squares (□) mark f{ai). 
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FIG. 2. Linear regression data for the evaluation of Eq. (Q) determining a(q) and f{q) for 
weakly coupled planes with 7 = 0.9, N = 48, W = 9 and q = -2(o), -!(□), 0(0), 1(A), 2(+). 
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FIG. 3. Wci'f) for weakly coupled planes as obtained from the "characteristic spectrum" (o) 
and from the system size dependence (A). The thin solid lines represent the two power law fits to 
the (o) data. The thick dashed line is the result of Ref. ^. The thick solid line is a combination of 
the isotropic resullH Wc = 16.3 and the perturbative exponent! /5 = 0.25 which fits the (A) data 
well. 
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FIG. 4. Wcij) for weakly coupled chains as obtained from the "characteristic spectrum" (o) 
and from the system size dependence (A). The thin solid line is a power law fit to the (o) data, 
the thick dashed line is the result of Ref. The thick solid line is the combination of the isotropic 
resulti Wc = 16.3 and the perturbative exponentl /? = 0.5. 



19 




FIG. 5. Singularity spectra for weakly coupled planes at 7 = 0.8 for the two system sizes 
AT = 18 (-) and iV = 42 (- -). The symbols distinguish W = 8(0), 10(A), 12(n) and indicate a{q) 
and f{q) for q = —2, —1,0, 1, 2 (from right to left). The error bars result from the linear regression 
(cp. Fig. ^) and the average over the different eigenstates (cp. Sec. IV A| ). The dotted line is the 
"characteristic spectrum" fc- 
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FIG. 6. N dependence of ao and ai for weakly coupled planes with 7 = 0.96 and 
W = 5.5(0), 6.5(n), 7.5(*), 8.5(0), 9.5(A), 10.5(+). 
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FIG. 7. Wave function close to the MIT for very weakly coupled planes with 7 = 0.99, N = 48 
and W = 4.5. Every site with probability larger than the average A^"^ is shown as a box 
with volume x A^. The 764 cubes with l^il"^ x N > VlOOO are plotted in white with black 
edges. The grey scale distinguishes between different slices of the system along the x-axis. The 
thick solid line is the logarithm of the summed probability amplitude for each plane perpendicular 
to the z-axis. Again only values above are shown. 
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FIG. 8. Density of states for weakly coupled chains (a) and planes (b) with 7 = 0.9 and N = 12. 
The data for = 5 and 15 agree with that of uncoupled chains or planes within the statistical 
fluctuations. The peak structure for = 1 is due to the small system size. 
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FIG. 10. Conductivity a for N = 12 and various anisotropies 7 and disorders W for weakly 
coupled chains (c), (d) and planes (a), (b) for the parallel (a), (c) and perpendicular (b), (d) 
direction. The perpendicular conductivity for 7 = 0.9 and = 15 is negligible while the parallel 
conductivity for 7 = 0.9 and W = I exceeds the range of the diagram (see Fig. P). 



25 



